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Abstract 

We present a mechanism for displaying the transmission property of the discrete 
Heisenberg ferromagnetic spin chain (DHF) via a geometric approach. By the aid of a 
discrete nonhnear Schrodinger-hke equation which is the discrete gauge equivalent to 
the DHF, we show that the determination of transmitting coefficients in the transmis- 
sion problem is always bistable. Thus a definite algorithm and general stochastic algo- 
rithms are presented. A new invariant periodic phenomenon of the non-transmitting 
behavior for the DHF, with a large probability, is revealed by an adoption of various 
stochastic algorithms. 



PACS numbers: 02.40.Ky; 05.45.Mt; 07.55.Db 



§1 Introduction 

There is current interest in displaying the properties of one-dimensional magnetic mod- 
els. The one-dimensional classical continuum Heisenberg models with different magnetic 
interactions have been settled as one of the interesting and attractive class of nonlin- 
ear dynamical equations exhibiting the complete integrability on many occasions ([l]-[7]). 
However, though the investigation of nonlinear spin chain systems is quite fascinating, 
few has been known in the case of more realistic physical lattice spin chains so far, es- 
pecially, for the following discrete (isotropic) Heisenberg ferromagnetic spin chain with 
nearest neighbor exchange interaction (DHF), 

Sn = S„ X (S„+i + S„_i), (1) 

where S„ = (s^, s^, s^) e TL^ with {s^)'^ + (s^)^ -|- (s^)^ = 1 and the dot stands for the 
time derivative. In fact Eq.([l]) comes from the Hamiltonian formalism: = {S,„, H} with 
the Hamiltonian function H = X]j=^oo Sj ■ Sj+i, where • denotes the inner product of 
vectors in R^, and the Poisson bracket {sj, s^} = 6jk Ylc=i ^abcSj, where 6jk is Kronecker's 
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symbol and Sabc is the 3-dimensional totally antisymmetric symbol. Physically, one would 
say the model ([T]) describes a system of classical spins and subjects to homogenous nearest- 
neighbor Heisenberg interaction. The standard continuous limit procedure performed on 
Eq.([T]) leads to the integrable Heisenberg ferromagnetic model: Sf = S x Sxx, which is an 
important equation on condensed matter physics (see, for example, [8]). Eq.([l]) is quite 
well known and there is little need to give here a detailed enumeration. It should be also 
emphasized that Eq.((T]) is widely believed to be not a completely integrable equation. To 
our best knowledge, there is no effective method in study of the dynamical behaviors of 
Eq.([T|) in literature. 

On the other hand, the geometric concept of gauge equivalence [U |T0] between inte- 
grable equations, which provides a useful tool in displaying solitonic dynamics, has been 
generalized to nonintegrable case in |1H [T2] . In [12] it is shown that the discrete nonlinear 
(nonintegrable) Schrodinger equation (AL-DNLS) is the discrete gauge equivalent to a 
nonintegrable discrete Heisenberg model and the transmission properties of the AL-DNLS 
equation ([I3]-[IH]) are completely preserved under the action of discrete gauge transfor- 
mations. In this paper, we would apply this geometric idea to transform the DHF ([1]) into 
a discrete nonlinear Schrodinger-like equation and utilize the nonlinear Schrodinger-like 
equation to study the transmission problem of the DHF ([T|) from mathematical point of 
view. We hope that this study will be helpful to reveal deeper quantum chaotic properties 
of the DHF ([T|) and useful in physical applications. 

This article is organized as follows. In section 2 we deduce a discrete nonlinear 
Schodinger-type equation which is the discrete gauge equivalent to the DHF in the cate- 
gory of the (discrete) Yang-Mills theory. In section 3 we display the transmission property 
for the DHF ([1]) with the aid of the discrete nonlinear Schodinger-type equation and an 
approximate linear stability analysis for the stationary version of the discrete nonlinear 
Schodinger-type equation is given for supporting the discussion of the transmission expo- 
sitions. Finally, in section 4, we close the paper with some conclusions and remarks. 

§2 Gauge equivalence 

Following the concept of discrete connection and associated discrete curvature introduced 
in [12], we would express Eq.([T]) geometrically as a discrete nonlinear equation with pre- 
scribed discrete curvature representation. We first convert Eq.([l]) into its matrix version: 




(2) 




4 + is: 



n 



,3 




Then we define a discrete connection by 



An = {Ln,Mn) 



(3) 
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where L„ = + M„ = i (l - ^^^j ^2i±|2i^ _ ^ z^-z ^ i+Sn^~iS^ ^-^j^ 

/ denoting the 2x2 identity matrix and z being a free spectral parameter. It is a 
direct and lengthy computation, in the similar way displayed in the appendix of |12j . 
that the corresponding discrete curvature {F^} is given by := L„ — M„+iL„ + 



-z''+z—z +z 



-(S„ • S„_i — S„ • S„+i)/. Here we have used the identity: Sn+iSn 
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-SnSn+i + (S„ • S„+i)2/ and similarly SnSn-i = -Sn-iSn + (S„ • S„_i)2/ in the com- 
putation. Thus we see that, if we set Kn = — — — (Sn ■ S„_i — S„ • S„+i)5„ + 
—2^+2-2 ^+z ^ (^s_^ . S„_i — S,„ • S„+i)I, Eq.([2]) is equivalent to holding the following pre- 
scribed discrete curvature representation: 

F^ ■.= Fn-Mn+lFn + FnMn = Kn, Vn. (4) 

Our aim in this section is to geometrically transform Eq.Q (resp., the connection ^) 
to a discrete nonlinear Schrodinger-like equation (resp., a new connection) by a discrete 
gauge sequence So the key point is to find such a For this purpose, for a 

matrix sequence {Sn\ solving Eq.([2]), we shall find {Gn\ to satisfy 

(T3 = —GnSnG~^ , Gn+i = i 9n | some (complex valued) qn, (5) 

\ -qn 1 / 

where ^-^ ~ q 1 ^ Pauli matrix. The first equation of ([5]) has a class of 

solutions of the form 

Gn = -==L==(CT3 - 5„)diag(F„, Fn), (6) 

v2(i - 4) 

where {-Fn} is free and to be specified latter. It is a straightforward verification that, in 
order for ([6]) to fulfill the second equation of ([5]) for some complex sequence {qn}, Fn is 
forced to satisfy the following iterated relation 



(l-4)(l-4+i)F„+iFn-i+(4-i4)(Sn+i+^4+i)^n+iF„~' = 2^(1 - si)(l - 4+i), Vn. 

So we may determine Fn progressively by this relation in n and hence have proved the 
existence of a desired gauge sequence {G„}. 

For the connection ([3]), we use {G„} to define a new connection {A^ = {L^{t, z),M^{t, z))} 

by L^{t,z) = Gn+iLnG-^ = Gn+i {^I + ^Sn)G-^ = I 1 and 



^ni^^ z) ~ GnGj^ ^ + GnMnG^ ^ — GnG^ ^ - 



-1 + 



z-^ 


zqn \ 






-z'^qn 


' J 


[l-Z^] 


1 911-1 


l+kn-l|^ 


il-z^] 







Here we have used the relations ([5]) in the computation. Since Ln, Mn satisfy the prescribed 
discrete curvature representation from Lemma 1 proved in [12] we know that L^, 
should fulfill 

L^-M^^,L^ + L^M^ = Gn+iKnG-\ Vn. (7) 
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It is a direct computation that 

1 Qn 



l.h.s. of dZD 



-qn 1 



U Z \ Z 



+ 1 



(z~z ^)qn . {z ^-z ^)q„ + (z-z ^)gn-l 



l+l9n|2 



l+|gn-iP 



8) 



where Cr, 



{Z ^-Z n-{z-Z l)|gn|^ (-Z ^+Z + ^)gngn-l J I 

^^-"'^lir^"')'^"^""-^ and 
l+|gn-ir 



(z-z^) + {z-z-^)\qn\^ , 



r.h.s. of (El) = i 



1 



1 



-1 



(9) 



Jn-1\ J y " + Z ")qn -Z + Z-" 

Here we have used the relation ([5]) and some identities displayed in the appendix of [12]), 
e.g., S„ • S„+i = ]''r!'^"L and S„ • S„_i = !~|'^"~^L , in the above computations. Thus, by 
substituting ([8]) and ^ into ([7|), the equation of the off-diagonal part of ([7|) leads to 



GuG-' = . _ . ^ , (10) 

where * are some expressions which cannot be carried out at this moment. On the other 
hand, at the same time we also have 

^nZ 



l.h.s. of dZD 



-qnZ ^ 

Gn+lG^\^ 



+ 



v-1 



z 



-1 



-1 



GnG„ 



-QnZ Z 
QnZ \ 



-QnZ ^ Z 



{z ^-Z^)qn _j_ {z-z^)qn + (z ^-z)gn- 



+i 



(z~z ^)qn . {z ^-Z ^)qri+{z~Z ^)gn- 



1+W 



+ 



LI) 



By substituting (fTTI) and dH) into the equation ([7]) and combining with (fTO]l . the diagonal 
part in this time implies 



-1 



■_2n _ „• Qn-l 



where a„ and satisfy 



/ q„qn- 



(12) 



(13) 



Obviously, the second equation of (|T3[) has a solution cr„ = q;„. By substituting ([12]) and 
([U]) into d?]), the off-diagonal part leads to {qn} satisfying the following discrete nonlinear 
Schrodinger-like equation 



qn + i 



qn+i 



n • qn , . qn—1 . , _ % „ 

^l- : TTT + l-r—^ + ("n - an+l)gn = 0, 



(14) 
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where solves the first equation of (fT3]) . This shows that the DHF ([T]) is the (discrete) 
gauge equivalent to the discrete nonlinear Schrodinger-like equation (fTl]l . 

We would point out that Eq.((TlD sets {^^ = {L^{t, z), {t, z))} as its discrete 
connection and has prescribed discrete curvature representation: = K^, where is 
given by right-hand-side of ([9]). Notice that the expression is to be zero at z = 1, it 
can also be proved conversely that the discrete nonlinear Schrodinger-like equation ()14p 
is (discrete) gauge equivalent to the DHF ([T]) by the gauge sequence {Gn} satisfying 
Gn+i = L^{t,l)Gn, Gn = M^(t,l)Gn, Vn (the existence of such a sequence {(?„} is 
because of the integrability of {^^}, or in other words its zero curvature representation, 
at z = 1). The details are omitted here and one may refer to [i2j for a reference. 

{an} appeared in Eq.([T3]) can be solved out from the first equation of as follows 
in different two approaches. One is obtained by iterating from initial data qq, g_i and oq- 

^^+^EU[TTlt:p-TTlt;^)Qk-m-2- n>0 
ao, n = (15) 



"l+knl 

The other is obtained by iterating from the boundary data at positive infinity: 

"n = i V ( — ^— 12 - -— -|^^ 12 1 Qk+iqk + i ^"^""^ Vn. (16) 

This reflects that the property of the nearest neighbor exchange interaction for Eq.([T]) col- 
lapses for its gauged equivalent equation (jl4p . since the n-th chain's exchange interaction 
in Eq. (ll4p relates to chains from —1-st chain to (n=b l)-th chain or all the chains with label 
^ n — 1. Though this implies that the dynamics of Eq. (ll4|) and hence the original Eq.([T]) 
will be complicated, the exposition of the approximate linear stability in the next section 
for the stationary version of Eq. (jl4p shows in interesting fact: the n-th chain term is related 
to its nearest neighbor exchange interaction in the stability analysis. This reflects that 
Eq. p4p does have the property of the nearest neighbor exchange interaction in some sense. 
We note that ( [H [71 |T2] ) the discrete nonlinear Schrodinger equations without nonconstant 
denominator terms (e.g. the AL equation: ig„-|-(g„+i-|-g„_i-2g„)-|-|g„p(g„+i-|-g„_i) = 0) 
are gauge equivalent to the (modified) discrete Heisenberg spin chain models with non- 
constant denominator terms (e.g., the Ishimori equation: Sn = 2 ^'^"s^^^^^'^s + 2]^^%-^). 
And, meanwhile, the above exposition indicates that the discrete Heisenberg spin chain 
model ([1]) without nonconstant denominator terms is gauge equivalent to the discrete 
nonlinear Schrodinger-type equation ()14|) with nonconstant denominator terms. This is 
an interesting duality phenomenon between discrete nonlinear Schrodinger-type equations 
and discrete Heisenberg spin chain models. 

Before we use ()14p to study the transmission problem of the DHF, let's give a general 
description of constructing solutions to the DHF ([T]), or equivalently (l2|), from those to 
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the discrete nonlinear Schrodinger-like equation ([HI). For a solution to Eq. p^ . let 
Go be a fundamental solution to 



,0 \, Go. (17) 

It is a direct verification that, by successive iteration, 

^ ^ f L^„i(t,l)G„_i, n>0 ^ f L^„i(t,l)...L^(t,l)Go, ^ > n 8) 

\ L^^i(t,l))-iG„+i, n<0 \ {L^^,{t,l)---L^{t,l)r^Go, n < 0^ 

solves G„+i = L^(t, 1)G„, G„ = M^(t, 1)G„, Vn. Therefore {S„ = -G^VsGn} is a solu- 
tion to the DHF ([2]) which corresponds to {qn} under the discrete gauge transformation. 



§3 Applications 

In this section, we shall use the gauged equivalent equation to study dynamical 
properties of the discrete Heisenberg ferromagnetic spin chain ([1]). We mainly focus on 
exploring the transmission property of the DHF ([T|) and its related linear stability analysis 
for the stationary version of the equation ()14p . 

A Transmission properties 

In this subsection, we study as a physical application whether the wave transmission 
property ([I3l[l6j) of the discrete nonlinear Schrodinger-like equation (fH|) and transfer 
it to that of the DHF ([T]) under the action of discrete gauge transformations. In order 
to do these, let's consider the following recurrence equation originated from the discrete 
nonlinear Schrodinger-like equation by setting qn{t) = exp(— ii?t): 



-^V'n + T— n Vi~^~r~\ 12 + ^i~n - i(an - a„+l)¥?n = O, (19) 



where 93^ is a complex amplitude independent of the time variable t, ii^ is a real parameter 
and a„ independent of the time variable t solves the first equation of (fT3j) . which reads 
from (fT6]l . 

m — 12 ~ m \2 "Pk+wk + 1——, [2 Vn. (20) 

We thus study the transmission problem of Eq. ()19p via a similar numerical way of the 
AL-DNLS lattice chain displayed in [T31[in] (please refer to their papers for details). First 
we note that Eq. (jl9p has a stationary solution as follows: 

ipn = Te*''", n G Z (21) 
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where the real parameter E satisfies the consistent condition 



2 cos K = 



2 + E{l + T^ 



) 



(22) 



We now consider the problem: A finite nonlinear segment O^n^A^ — lof length in 
the nonnegative stationary regime is embedded in a nonlinear chain {ipn}nez satisfying 
(dHD with 



where (j22p is fulfilled. This can be regarded as that an incident plane wave Roe^'^ on the 
left (n = 0) induces a reflected plane wave Re~^'' on the left and a transmitted plane 
wave Te***" on the right [n ^ A^). Rq is called the amplitude of the incoming wave, R the 
amplitude of the reflected waves and T the transmitted amplitude at the right end of the 
nonlinear chain; k is called the out-coming wave number; |-RoP ^-nd |Tp are also called the 
in-coming and the transmitted intensity respectively. The medium is completely nonlinear, 
thus the transmission coefficient T as a function of Rq may not uniquely determined. If 
this is true, according to the sense made in |17[|13j . it is called bistability. For the forward 
transmission problem of the AL-DNLS lattice chain, there occurs exactly the bistability 
phenomenon ( ]1?>\ ITB] ) , which leads Delyon et al [13] to consider the backward transmitted 
problem for the DNLS. They proved that the pair (k, |r|) initializes the incident intensity 
|i?o| completely and then displayed the transmission behavior of the DNLS in the (k, |r|) 
plane. 

Now we consider the similar backward transmission problem for solutions with the 
type ([^5]) to the present equation Following Wan and Soukoulis [T3], we use polar 

coordinates for that is = -R„e*^". For a given pair (k, |r|) (without loss of the gen- 
erality, we may assume T > in this paper), we try to iterate Eq. ()19p from n = A^ to — oo 
successively and to determine the amplitudes {Rn-i, • ■ ■ , Ro) and phases {6n-i, • ■ ■ , ^o)- 
The existence of solutions with the type (|23|) to Eq. (|19|) will be shown below in the algo- 
rithms. In order to further support the existence of such solutions, an approximate linear 
stability analysis around the the stationary solution (I2ip will be given in the next section, 
which further reveals the linear stability behaviors of the equation ()19p . In this iterating 
process, we find that there are two choices in determining the amplitudes {Rn-i, • • • , -Ro) 
in each step. This implies that, not like the transmission problem of the AL-DNLS, the 
pair (k, T) does not initialize the incident intensity \Rq\ completely, or in other words, the 
backward transmission problem for the equation ()19p with solutions of type ()23p is still 
bistable. 

We first introduce a definite algorithm such that the pair (k,T) initializes the am- 
plitudes (i?Af_i, • • • , Rq) and phases {On-i, ■ ■ ■ ,9i) completely. In fact, such an algorithm 
is designed as follows. From ipn = Te*''", n N we get a„ = i ^^[""^""^p , n ^ N from 




(23) 
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(j20j) and the equation p9j) for n ^ + 1 is now equivalent to (j22|) . The equation (|T9|) for 
n = is rewritten as 

Tsin(K) - -^±i^fi^_i sin(A^^) = 0, 

where A0„ = On — On-i (we remark that ^tv = i^^), which leads Rn-i to satisfying 
T = y:^^? — and ^at-^i = k(A^ — 1). Since the quadratic equation T = j^rx in 
X just has two self-reciprocal solutions: x = T and x = we then determine Rn-i 
according to the rule: to choose Rn-i > 1 when < T < 1 and Rn-i < 1 when T > 1 
in order to avoid getting the stationary solution (f2T|) . Next, from (f20l) we see that, for 

n<N,an = i Ek=n+i (iq^TTF ~ TT^) '^kVk-i + h+I^T-iV ' '^^^'^ ^^"^ equation ((121) 
for n < is is equivalent to 



(24) 



^"-^ rn.fA^' 1 - f g(l+4)+2 i?n+i( j;^-fi|+2R„+2R„i?^^,) 
T+R^ cos(A(y„ j - ^ I cos(A6y„+i J 

^iTfijs- - iri^) cos(A6'fc+i)|, 
-^^sm{A9n) = y|^cos(A0„+i). 

The above exposition indicates that Rn-i satisfies 
R 

= ^/fr.h.d of the first eq.(l2iD)2 + fr.h.d of the second eq.ipill)^" 

and thus also has two self-reciprocal solutions. Our rule is that in the second step and 
the following steps, we always choose the root with Rn > 1 for the quadratic equations. 
According to this algorithm, we may uniquely determine the amplitudes and phases from 
n = A^ — Ito— oo successively and especially the incident intensity |iio|^- This procedure 
also indicates the existence of solutions of type ([^5]) to the stationary equation . 

If the resulting incoming wave intensity |i?oP is of the same order of (the reciprocal 
transmitted intensity |Tp) independent of when < T < 1, or is of the same order 
of |rp when T > 1, we say that the nonlinear chain with wave number k and outgoing 
intensity |Tp is to be transmitting. If Rq appears to be a rapidly increasing function of 
A^, we say that this nonlinear chain is to be non-transmitting. Figs. 1(a) and (b) below 
display the transmission behaviors in the {k,T) parameter plane with the chain length 
A^ = 110 and 150, respectively, representing region of transmitting (white) and non- 
transmitting (black) behaviors. In particular, Fig.l (a) and (b) numerically show the 
whole-scope of the transmitting behaviors for k G (1,3.5) and T G (0.5,2). The non- 
transmitting regions form some strange patterns. The enlargement of one of these patterns, 
(k, r) E (1.3, 1.7) X (1, 1.4), is further shown by Fig. 1(c). It is seen from the figure that 
the twisting boundaries between the transmitting region and the non-transmitting region 
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in the parameters plane always display spiral, cantor-like, and fractal structure, which is 
very analogous to the structure of the standard chaotic attractor. Furthermore, from the 
Fig.l (a) and (b), we see that the transmitting behaviors only depend seriously on the 
parameters k and T but do not depend sensitively on the chain length N . 
Figure 1 around here. Please find it at the end of the paper. 

Besides the definite algorithm, we come to consider a stochastic algorithm. In such 
a stochastic algorithm, the determination of the bi-valued Rn-i should be the same as 
that in the definite algorithm in order to avoid getting a stationary solution; however, the 
latter determinations of i?„ (0 ^ n ^ — 2) in each step are based on the generation of a 
random number ^, i.e., Ri = R- ^ c; Ri = R.^ otherwise. Naturally, our first choice of 
the generation of the random number ^ is the normal distribution and the value c is taken 
as zero, the expectation of this distribution. It could be seen from Fig. 2 (a) that non- 
transmitting regions (black) periodically appear with respect to n, here the chain length 
is still to be = 150. Also, this kind of phenomena could be somewhat found in Fig.l 
(a) and (b). Secondly, a Poisson distribution with an expectation 5 is adopted to produce 
the random number ^ and c is taken as 5. In this case, a similar periodic non-transmitting 
behavior is displayed in Fig2.(b). Actually, thousands of numerical simulations verify a 
fact that no matter what kind of stochastic distribution and chain length are taken, the 
non-transmitting behavior is so prevalent around the values kj. ~ 1.57+ ^ {k = 0,1, ■ ■ ■). 
This periodic phenomena could be regarded as an invariant non-transmitting behavior for 
Eq.dlS]). 

Figure 2 around here. Please find it at the end of the paper. 

The transmission problem associated with the discrete Heisenberg ferromagnetic spin 
chain ([1]) (or equivalently (l2|)) is now proposed as follows. There is a semi-infinite nonlinear 
matrix wave chain {Sn = —G^^a^Gn} embedded in a nonlinear chain of the DHF ([2]), 
where Gn looks like (from (fT8|l ): 

\ -qn-i 1 y \ -90 1 / 

with Go satisfying (fT7|) and qn = i/^ne"*^*, where ipn is independent of t and given by 
()23p . This nonlinear matrix wave chain is from the left towards the right, where they are 
scattered into reflected and transmitted parts. Similarly, Ro,R present the amplitudes of 
the incoming and reflected waves and T the transmitted amplitude at the right end of 
the nonlinear chain under consideration. |-RoP and jTp are also called the incoming wave 
intensity and the transmitted intensity of the nonlinear matrix chain respectively. As we 
have proved in the previous section that the DHF ([T]) (i.e. the DHF ([2])) is the discrete 
gauge equivalent to the discrete nonlinear Schrodinger-like equation (fH]) and verse visa, 
we see that the transmission properties of the equation (fT9]) just displayed, i.e., the above 
definite algorithm and stochastic algorithms, are transferred completely to the trans- 
mission problem (I25p of the DHF ([2]) under the action of discrete gauge transformations. 
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Thus the nonhnear matrix wave chain (j25p of the DHF ([2]) has the same transmission 
behaviors as those of nonhnear chain l\23h of Eq. ()19p . Fig.l(a),(b) and (c) are completely 
suitable in displaying the transmission behaviors of the finite nonlinear matrix chain of 
the DHF ([2]) in the parameters {k,T) plane with chain length N = 110 and 150 under 
the definite algorithm and so do Fig. 2(a) and (b) under stochastic algorithms. Further- 
more, from the Fig.l (a) and (b), we see that the transmitting behaviors for solutions 
{Sn = —G~^asGn} with {Gn} satisfying I^E\i to the DHF ((T|) only depend seriously on 
the parameters k and T but do not depend sensitively on the chain length A^. This shows 
that the important phenomenon of bistability occurs still in spin magnetic physics. The 
periodic phenomenon of the non-transmitting regions displayed in Fig. 2 provides a new 
and interesting transmission property of the DHF ([T]) . 

B Stability analysis 

In order to support the exposition in subsection A and show the existence of stable so- 
lutions to (fT9]) around the stationary solution (f2T]) . we shall give an approximate linear 
stability analysis for Eq.p9|). In what follows, we focus on the locally linear stability of the 
stationary solution (|2ip in the reverse iterative procedure through the linear variational 
method [IHIEO]. For this purpose, we naturally import the formulas as follows: 

Rn=T + Tn, 9n = nK + 5n, n = N, N - 1, ■ ■ ■ ,2,1, 

where quantities t„ and 5„, respectively, are regarded as small perturbations of the mod- 
ulus and argument of the stationary solution (j2ip . Substitute the above formulas into 
Eq. (|24p and then expand those nonlinear terms in the vicinity of T and uk. A tedious 
calculation thus gives the following equations at 0(rj) and 0{6i) {i = n — l,n,n + 1): 

' (l-r^)cOSK 



lyj^i i^n ~ ^n-l) 
(l-T^) siriK 
(1+T2)2 '^n-l 

TCOSK f X X \ 

where n = A^, A^ — I,-- - ,2,1. Clearly, combining the second and the fourth equations in 
([26D leads to 

^n+l — (^n = 0, 

which shows an identity of the series {6n} at 0(5„). Hence, the linear stability of the 
argument of the stationary solution in the reverse iteration procedure belongs to the stable 
case. However, the nonlinear stability of this argument becomes a critical case, depending 



1-T2 

(1+T2)2 

I T ^ 
I T ^ 



cos K 



T2 cos k 
1+T2 



I COSK _ 2T(1-T2)cosk " 
T{l+T2) (i+r2)2 

(1~T2)cosk 2rcosK I 2T(1-T2)cosk " 

^''^^^^^''^'^ 1+T2 (1+T2)2 



4T 

'(1+T2)2 



T(l+T2)2 

COS K sr^N—l { 
+ (1+T^)2 Luk=n+\Vk - Tfc+l 

(l-T2)sinK 
T cos K I X X \ 



(26) 
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on the further calculation at 0{df). Due to the tediousness of notations, we omit the 
calculation here. Analogously, combining the first and the third equations in ([26]) yields: 



Tn-l 



1 + 



2T' 



1-T4 (l-T4)cOSK 

Hence, it follows that for n = A^, • • • ,1, 



2r 



1 -T^ 



TN- 



Tn-1 



2T' 



-I N-n+l 



1 - 

2T^ 



(1 

N-n 

E 



j=Q 



T^) cos K 

1 + 



2r2 



1 - r4 (1 - cos K 



'Tn — XN,nTN- 



Consequently, the linear stability of the modulus of the stationary solution is stable in the 
reverse iteration provided that \N,n is uniformly bounded for all N and n. More precisely, 
the reverse stability holds when the following inequalities are satisfied: 

2T2 



1 + 



1-T4 (1-T4)cosk; 



< 1. 



(27) 



In addition, notice that r„ = A^qTo. This implies that the linear stability of the normal 
iteration procedure is stable when A„.o ^ for all n. Thus, this linear stability is valid 
only provided 

fi^^ ik = 0,±l,±2,---). 

Obviously, this necessary stability condition is also included in the condition ()27p . As 
shown in Fig. 3, parameters selected from the black regions violate the inequality ()27p . 
Those diamond-like unstable regions periodically appears in a period 2tt with respect to 
the parameter k. 

Figure 3 around here. Please find it at the end of the paper. 

The relation of the linear stability analysis in this subsection and the transmission 
exposition in the subsection A is that the non-transmitting regions in Figure 1 or Figure 2 
are always contained in the unstable regions in Figure 3 (b). This also indicates that the 
transmission exposition for solutions with the type p3p to Eq. (ll9p is a more accurate de- 
scription of the linear stability analysis for Eq. ()19p , which exhibits some chaotic dynamics 
of the Eq.([Tl]) and hence those of the DHF 



§4 Conclusion 

We present a mechanism for displaying the transmission property of the discrete Heisen- 
berg ferromagnetic spin chain via a geometric approach. More precisely, we first transform 
the DHF (pP) into the discrete nonlinear Schrodinger-like equation (|14p with the aid of the 
(discrete) Yang-Mills theory. Then we use the stationary version of the discrete nonlinear 
Schrodinger-like equation ()14|) to study the transmission problem of the DHF ([1]). In this 
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procedure, we show that the determination of transmitting coefficients in the backward 
transmission problem is always bistable. Thus a definite algorithm and general stochastic 
algorithms are presented. The corresponding transmitting behaviors in the parameters 
(k, r) plane are shown by Figure 1 and 2. A new invariant periodic phenomenon of the 
non-transmitting behavior for the DHF, with a large probability, is revealed by an adop- 
tion of various stochastic algorithms. By the way, an approximate linear stability analysis 
for the stationary version of Eq. (|19|) is also given for supporting the discussion of the 
transmission expositions. We remark that though the gauged equivalent equation (jl4p 
with (I13p looks much complicated and we cannot give a direct physical application at 
present stage from it, it is an effective way to utilize it to study some dynamical properties 
of the DHF ([T]) . We believe that the discrete nonlinear Schrodinger-like equation p4|) can 
help us to reveal much more and deeper quantum chaotic properties of the DHF ([1]). 
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Fig. 1: The transmitting behavior in the {k,T) parameters plane {k,T) e 
(1,3.5) X (0.5,2) (a) with N = 110 and (b) with iV = 150, and in an en- 
larged region (k, T) e (1.3, 1.7) x (1, 1.4) (c) with N = 150. The boundaries of 
transmitting (white) and non-transmitting (black) regions exhibit fractal-like 
structure. . 
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Fig. 2. The periodic phenomenon of the non-transmitting behavior in {k, T) plane 
by adopting of the normal distribution (a) and a Poisson distribution (b). Here the 
chain length N = 150. 
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Fig. 3. The linear stability parameters regions (white for stable region and black 
for unstable region) when (k,T) is in (0, 13) x (0,2) (a) and in (1,3.5) x (0,2) (b), 
respectively. 
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